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First of all, we mention the situation of empirical analyses on 3rd order BEC (Bose-Einstein 
Correlation) at RHIC. Second, we introduce several theoretical formulae / approaches. Third we 
present our analyses of data in Au+Au at 130 GeV by STAR and preliminary data in Au+Au at 
200 GeV by PHENIX Collaborations. Our results also contain analyses by means of core-halo 
model. Finally, we estimate that the volume of interaction in Au + Au collisions at 130 GeV is 
500 fm^, which is compared with V — ^long^out^side ^ 300 fm^^ in Pb + Pb collision at 2.76 TeV 
by ALICE Collaboration. Moreover, usefullness of empirical analyses on {27t^)7t^ and (27r^)7r+ 
combinations at RHIC and LHC energies is remarked. 
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1. Situation of Empirical analyses on 3rd order BEC in Au+Au collision at RHIC 

As shown in Table [l], STAR and PHENIX Coll. have reported their analyses at 130 GeV and 
200 GeV, respectively. In their analyses, STAR Coll. [|I|] used all combinations of three momentum- 
transfers ^J^j < (Inside of a globe in Fig. |l|), 

2L,3 ='?12 + '?23+^31• 
On the other hand, PHENIX Coll. used data on diagonal line of a cube in Fig. [Tj, 

q3 = {qn) = {q23) = (qsi) i-e. eL,3 = 3^3- 
Of course this relation holds, as the number of data increases. (• • •) is an average value. 

Table 1: Situation of analyses on the 3rd order BEC. Notice two empty columns. 
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(1.2) 





STAR 


PHENIX 


130 GeV 


raw and corrected data by Qi,n, 




200 GeV 




preliminary raw and corrected data by ^3 





yqf2 ^^'^i'2 V3q3 

Figure 1: Data ensembles of STAR (left) and PHENIX Coll (right). 

Here we compare two kinds of data. In Figs. ^ data by PHENIX Coll. are rearranged by 
V^q3. Coincidence among data by STAR and PHENIX Coll is fairly good. Error bars in raw data 
by PHENIX Coll are smaller than those of STAR Coll. 

2. Several theoretical formulae 

In many analyses on BEC, the following formulae based on plane wave function are used. 



i+A£^,-(«e,v)^ + 2Ai-5^.-°-5(«e3)^ 



(2.1) 
(2.2) 



In laser optical (LO/GL) approach, the following formulae with a degree of chaoticity p have 
been proposed [0] and utilized. 



(2.3) 
(2.4) 



where Ej^ = exp(— (Gaussian form) and/or = exp(— /?y^2^) (exponential form), and 
E^g = exp{—R^Q^) and so on. 
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Figure 2: Comparisons of data by STAR and PHENIX Coll. 

Third we explain the formulae by Coulomb wave function including the degree of coherence 
X and the interaction range 7? [|, |, |, ^. The two-body Coulomb wave function is well known as, 



^iflk^{xiXj)=T{l + i^ij)e'''^''l^e'''•y''^^F[-i^iJ, 1; /(^^^r^y - • r^y)], 



(2.5) 



where, r^y =Xi —xj, kjj = {kj — kj)/!, rjj = |r,y|, kij = \kij\ and T],y = eiejUij /kjj. jXij : reduced mass 
of m; and my, F[a, b; x] : confluent hypergeometric function, r{x) : Gamma function. 

Using of Eq. (2.5), the 2nd order BEC with A and Gaussian form for p(x,) is calculated as, 

!=I 

=n/ pi^M^^i \ {\^,kS^^^ -^2)1^ +1^^^^^ (^2, -^i)i^) 

+ARe(^i/A[^^^(xi, X2)i//[*^^(x2, xi)) . , (2.6) 
The 3rd order BEC is computed based on the 3-body Coulomb wave function in ||^, 

"^f = '/^&,(^l'^2)V^g^^(x2,X3)l/^f;^/x3,Xi), (2.7) 

Hereafter, we use the following expression; y/^'^ (x/,xy) = e'^^/^'^'^'''^^^ (r,y). Notice that the 
numerical factor "2/3" in the exponential function is important [^, ^. 



\t\lp{x{)d\,\Y^A{j)\\ 

,=i J y=i 



where 



A{2) = A23 = wlkM, ^sHlkM^ ^2)'/^r. {X2, xi) 



(2.8) 



(2.9) 
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Figure 3: Diagrams of the 3rd order BEC for A(l) ~ A(6). 



Aijk is reflecting the permutations of particles, /, j, k in Fig. |3|: Therein A(3) = An, A(4) = A123, 
A(5) = A132 and A(6) = A13. In the plane wave approx., we have the correct expression, where 
A(3) ~ A(6) are skipped [|, ^, 



PW 



A(l) = Ai — 



J{2/?,){kn-rn+k2i-r2i+hi-rii) _ i(ki-Xi+krX2+kyXi) 



A(2) = A 



23 



j]W i(2/3)(/ci2-ri3+fe3-r32+^3i-f2i) _ j(ki-x^+krXi+krX2) 



(2.10) 



Combining A(l) ~ A(6) in Fig. |3|, we obtain F\ as: (The other formulae F12 ~ ^123 and F132 are 
given in [|, ^.) 



1 



Fi = -[AiA*+Ai2At2+A23A^3+Ai3A*3+Ai23At23+Ai32At32]. 



(2.11) 



cfl / p{xi)d\i [Fi + 3AFi2 + 2A^Re(fi23] 



C 



(2^/371/?: 



2\3 



2/?2 12^1 + 3^2 



Fi + 3AFi2 + 2A5Re(Fi23: 



.(2.12) 



where = X2 — xi, = ^3 ~ (wii^i +m2X2)/M2, (^3 = (mixi +^2X2 + ?n3X3)/M, M2 = mi +/«2 
and M = mi+m2 + m^. 

Finally, it is known that the core-halo model is a useful model. However, explicit expressions 
are skipped here. See our studies in [|, 1,11,0]. See also [§]. 

3. Analyses of data by STAR Coll and PHENIX Coll 



3-1) For Coulomb corrected data, we employ the conventional formulae, Eqs. ( |2. 1] ) and ([2.2|). 
Ours are given in Table ^ and Fig. ^ It is interesting that R2n ~ Rj,n ~ 8.5 fm for data by STAR 
Coll. 



Table 2: Analyses of corrected data by STAR and PPHENIX Coll (Eqs. (2.1) and (0.2).) 
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Z^/n.d.f. 


R [fm] 
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271" 


8.75±0.31 


0.58±0.02 


23.0/25 


4.77±0.04 


0.39±0.01 


178/40 


371- 


8.26±0.39 


0.50±0.02 


1.88/35 


6.92±0.82 


0.34±0.10 


6.5/14 



3-2) Corrected data with ^3 > 0.02 GeV by PHENIX Coll are analyzed by conventional for- 
mula / LO approach. Our results are shown in Figs. ^ and ^, and Tables ^ and ||. 

3-3) Raw data with ^3 > 0.02 GeV are analyzed by the formulae of Coulomb wave function 
(Eqs. (2.6) and (2.12)). Our results are also shown in Fig. ^ and lower-part of Table ||. 
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Figure 4: Analyses of corrected data by STAR Coll. and PHENIX Coll. Eqs. (gjj) and (gj) are used. 
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Figure 5: Corrected data are analyzed by Eq. (2.3) and (2.4) 
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Figure 6: Analyses of raw data by PHENIX Coll. Eqs. (2.6) and (2.12) are used. 



3- 4) Using two formulae in the core-halo approach (with Gaussian source function, the fraction 
of core part /c and the degree of coherent pc in [Q, ^, ^) we obtain Fig. ^ In raw data, there is 
no over-lapping region. On the other hand, in corrected data, we observe very narrow over-lapping 
region. The reason of the wide region of 3?! BEC is due to large error bars of corrected data 

4. Summary 

4- 1) From raw data as well as Coulomb corrected data in Au-i-Au at 130 GeV by STAR Coll., 
we get the following interaction ranges R2n = 8-7 fm and Rt,]i = 8.3 fm, and can estimate 

V =Rl^^5mfm^. (4.1) 

This value is compared with that of ALICE Coll [^], V = /?iong^out^side ~ 300 fm^ at dNch/drj = 
1500 and kr ~ 0.3 GeV in Pb-^Pb at 2.76 TeV 

4-2) On the contrary, from corrected data at 200 GeV by PHENIX Coll., we obtain the ranges. 



by utihzing Eqs. (2.1 ) and (2.2), = 4.8 fm (A = 0.39) and Rsj, = 6.9 fm (A = 0.34). From raw 
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T able 3: Analyses of data by PHENIX Coll. Eqs. (|I6|) and (p^ are use d 
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Figure 7: Analyses of data by core-halo model. 

data, we have smaller interaction ranges = 3.8 fm and Ri,ji = 5.8 fm. The interaction ranges of 
data at 200 GeV by PHENIX Coll. are smaller than those of STAR Coll. At present, it is difficult 
to draw concrete physical picture for Au-i-Au collision at 200 GeV. Then we are waiting for final 
empirical analyses by PHENIX Coll. 

4-3) Moreover, we also eager for empirical analyses of (27r+)7r^ and (271^)71+ combinations 
at RHIC and LHC energies. See [0, [Tl|]. 
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